Evidence for incommensurate spin fluctuations in Sr2Ru04. 
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We report first inelastic neutron scattering measurements in tlie normal state of Sr2Ru04 that reveal 
the existence of incommensurate magnetic spin fluctuations located at qo = (±0.67r/a,±0.67r/a,0). 
This finding confirms recent band structure calculations that have predicted incommensurate mag- 
netic responses related to dynamical nesting properties of its Fermi surface. 



Being the only example of a non-cuprate layered per- 
ovskite superconductor [||, Sr2Ru04 has attracted con- 
siderable attention despite its rather low critical temper- 
ature , Tc '^l K Its normal state is characterized as 
an essentially two-dimensional Fermi liquid and the co- 
herent interlayer transport settles in at low temperature 
only 1^. The susceptibility is most likely dominated by 
an enhanced Pauli spin susceptibility, x- Meanwhile, the 
Sommerfeld coefficient, 7, in the specific heat is enhanced 
by a factor 3.5 with respect to band structure calculations 
||,|. It yields a Wilson ratio, Rw (~ x/7), of 1.7. This 
value indicates that the enhancements in both suscepti- 
bility and electronic specific heat can be ascribed to the 
same origin: most likely correlations among electrons . 

Noticing that SrRuOa is ferromagnetic (FM), it has 
been conjectured that Sr2Ru04 is close to a FM instabil- 
ity as well m. This assertion is corroborated by micro- 
scopic calculation of magnetic properties of ruthenates 
Since FM fluctuations disfavor both s- and d-wave 
superconductivity, it has been suggested that supercon- 
ductivity in Sr2Ru04 should possess p-wave symmetry 
(triplet pairing) . Conventional local density approx- 
imation (LDA) calculations give a correct Fermi- 
surface topography, probed by de Haas- van Alphen mea- 
surements as well as the magnetic enhancement due 
to Stoner exchange enhancement, although the mass 
renormalization cannot be explained within LDA calcu- 
lations. In the superconducting state, the ^°^Ru nuclear 
spin lattice relaxation rate, 1/Ti exhibits a sharp de- 
crease without a coherence peak (Hebel-Slichter peak) 
just above Tc, supporting the idea that an anisotropic 
pairing is effectively realized in Sr2Ru04 [|ll|. In ad- 
dition, the spontaneous appearance of an internal mag- 
netic field below the transition temperature, reported by 
muon spin rotation measurements (fiSR) [ p^ , and the 
absence of ^'^O Knight shift modifications below Tc 
point towards the triplet p-wave superconductivity. How- 
ever, few experiments have really probed the exact nature 
of the spin fiuctuations. Only the observation of a simi- 
lar temperature dependence for ^°^Ru 1/TiT and for ^''O 
1/TiT in the NMR experiments by Imai et al. has 
suggested that spin fluctuations are predominantly FM 
in origin. 

The determination of the antiferromagnetic order in 



the closely related compound Ca2Ru04 [|l^,|T^ has sug- 
gested that the picture of a near-by FM instability 
in Sr2Ru04 is too simple. Furthermore, recent calcu- 
lations which take into account the particular topol- 
ogy of the Fermi-surface, have predicted a sizeable 
magnetic response at the incommensurate wave-vector 
(27r/3a,27r/3a,0) |10|, i.e. far away from the zone- 
center. The enhanced susceptibility arises from pro- 
nounced nesting properties of the almost one-dimensional 
<i-xz,yz bands. Mazin and Singh discuss the possibility of 
a competition between p-wave and d-wave superconduc- 
tivity in Sr2Ru04 @- 

In this letter, we report first inelastic neutron scatter- 
ing (INS) measurements performed on single crystals of 
Sr2Ru04 in the normal state. Our data reveal dominant 
magnetic scattering at the incommensurate wave vectors 
qo=(±0.67r/a,±0.67r/a,0), i.e. very close to the positions 
predicted by the band-structure calculations. The rele- 
vance of these findings for the mechanism of supercon- 
ductivity in Sr2Ru04 will be discussed. 

Most of the INS measurements presented here have 
been carried out on a single crystal of cylindrical shape 
(4mm in diameter and 35mm long) grown by a fioat- 
ing zone method. The sample exhibits the supercon- 
ducting transition at Tc ~0.62 K. The single crystal was 
mounted in an aluminum can and attached to the cold 
finger of a closed cycle helium refrigerator. The INS ex- 
periments were performed on the triple axis spectrom- 
eters 2T (thermal beam) and 4F2 (cold beam) at the 
Laboratoire Leon Brillouin, Saclay, France. These spec- 
trometers use neutron optics that focus the beam to the 
sample, with a resulting gain of neutron flux that proved 
to be crucial for these experiments. The experimental set 
up incorporates PG002 monochromator and analyzer and 
14.7 meV fixed final energy. A pyrolytic graphite filter 
was inserted into the scattered beam in order to remove 
higher order contaminations. Data were taken within 
the scattering plane spanned by (1,0,0) and (0,1,0) di- 
rections. Some additional measurements were performed 
using several smaller single crystals with higher transition 
temperatures, rc=1.4-1.5 K; these experiments have re- 
vealed similar signals. Throughout this article, the wave 
vector Q=(H,K,L) is indexed in units of the reciprocal 
tetragonal lattice vectors 27r/a — 27r/6 — 1.63 A^^ and 
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FIG. 1. Constant-oj scans performed at huj—6.2 meV 
around Q=(1.3,0.3,0) along the (0,1,0) direction: T=10.4 K 
(•), T=295 K (o). 

2tt/c = 0.49 A"^ (I4/mmm space group) 0. 

Figure |l| shows representative constant-u; scans taken 
in the (H,K,0)-plane: at fiuj=6.2 meV and around Qo 
= (1.3,0.3,0) along the (0,1,0) direction. The scan at 
10.4 K shows a sharp maximum of intensity peaked at 
Qo=(1.3, 0.3,0) on top of a smooth background. At room 
temperature, this sharp peak has almost disappeared. 
The horizontal bar indicates the spectrometer resolution. 

At 10.4 K, several constant-w scans, with 6.2 meV 
energy transfer and performed along different direc- 
tions ((1,0,0), (0,1,0), (1,1,0), (1,-1,0)) have revealed the 
existence of comparable peaks at Qo =qo-l-G, where 
qo=(±0.3,±0.3,0) = (±0.67r/a,±0.67r/a,0) and G is a 
zone-center or a Z-point (001) in the (HKO)-plane. The 
fit of the data to a Gaussian profile incorporating exper- 
imental resolution function demonstrates that the peak 
intensity is isotropic with an intrinsic q- width (FWHM), 
Ag=0.13 ± 0.02 A^^. 

The interpretation of the scattering at qo as magnetic 
in origin is supported by the large number of points in 
reciprocal space where it has been observed. Further, the 
lowest phonon frequencies at qo are above 12 meV [^T). 
In addition, in contrast to a phonon-related scattering 
that increases at large \Q\ or with temperature, the scat- 
tering at qo decreases both at large wave vector (Fig. ||) 
and at high temperature (Fig. [|). These different points 
establish the magnetic origin of the scattering observed 
around qo. In contrast, in spite of several attempts, no 
sizable FM spin fluctuations have been observed. 

In a paramagnetic state, the magnetic neutron cross 
section per formula unit can be written in terms of 
the imaginary part of the dynamical spin susceptibility, 
X"(Q,C^), as [^9|j20); 
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FIG. 2. Magnetic intensity, measured at T=10.4 K and 
hu}=6.2 meV as a function of \Q\. For each point, the corre- 
sponding wave vector, (H,K,L), is also reported. The full line 
corresponds to the square of the Ru'^ magnetic form factor. 



tor and g ~2 is the Lande factor. The intensity of 
the scattering can be reasonably well described by the 
squared magnetic form factor of the Ru+-ion ||l8| (note 
that the magnetic form-factor of Ru''"'" is not available) 
after correction for geometrical factors related to the un- 
favorable shape of the sample, see Fig. 2. According to 
our measurements, the q-dependence of x" is given by: 
X" (Q, u;) = x" (90, CO) exp[-4 ln(2)(Q - Qo) VAg^]. 

The Fermi surface in Sr2Ru04 is formed by three 
sheets |l^: one, related to the 4(ia;j,-orbitals is quasi-2D, 
whereas, the two others, related to 4:dxz.yz orbitals are 
quasi-lD. The ID-sheets can be schematically described 
by parallel planes separated by g=±27r/3a, running both 
in the x and in the y directions. These peculiarities 
give rise to dynamical nesting effects at the wave vec- 
tors k={q,ky), k={kx,q) and in particular at q=(g, 9). 
The nesting effects become dominant when calculating 
the bare spin susceptibility of a non interacting metal 
|l0| , given by the Lindhard- function [IiqI : 



Xo(g,t^) = -2^1^ 
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fk 
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where rg =0.292 barn, f (Q) is the magnetic form fac- 



where e -^0, fk is the Fermi distribution function and 
£k the quasiparticle dispersion relation. Our INS are 
in very good agreement with the predicted four spots 
of magnetic scattering situated at q=(±27r/3a,±27r/3a) 
]lOt . In the experiment the incommensurate mag- 
netic responses are actually observed slightly away, at 
qo// = (±0.67r/a,±0.67r/a), which is most likely related to 
details of the band-structure jlOj . 

Let us now consider the energy dependence and mag- 
nitude of x"(9oj'^)- At T=10.4 K, constant-w scans 
have been measured at Q=(1.3,0.3,0) along the (0,1,0)- 
direction for different transferred energies between 2.4 
and 12 meV. The magnetic response always displays 
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FIG. 3. Energy dependence of the imaginary part of the dy- 
namical magnetic susceptibihty at Qo=(1.3,0.3,0) as obtained 
from energy scans (o) and constant energy scans around Qo 
along the (0,1,0) direction (•) (see text). 

a Gaussian profile, located at qo with an energy in- 
dependent q-width, on top of a constant background. 
In addition, two energy scans have been performed at 
Q=(l. 3,0.3,0) and at Q=(l. 3,0.46,0), the latter provid- 
ing a background reference. These measurements allow 
us to determine the energy dependence of the magnetic 
response at qo from 1.5 to 12 meV. The analysis could not 
be extended to higher and lower energies due to the con- 
taminations by phonon ||l7|| and elastic incoherent scat- 
tering respectively. Using Eq. ([^) , the magnetic intensity 
has been converted to the dynamical spin susceptibil- 
ity x" after correction by the thermal population factor 
and the squared magnetic form factor reported in figure 
^. We have then calibrated x" in absolute units against 
acoustic phonons, according to a standard procedure . 
X"(w, Qo), whose energy dependence is reported in abso- 
lute units in Fig.||, slightly increases up to 7 meV and 
then almost saturates. This energy dependence can be 
parameterized following linear response theory: 



x"{qo,uj) = x{qo,0)- 



Tlo 



Y2 



(3) 



where F is a damping energy of 9 meV and x'(90;0) = 
180 ^^.eV~^ corresponds to the static spin susceptibil- 
ity at qo. It is worth emphasizing that x'(90jO) is 6 
times larger than that at Q=0, i.e. the uniform sus- 
ceptibility: X = x'iQ 0,0) = 30 ^I%.cy-^ (~ lo-^ 

emu/mole) |l|-^. In La1.86Sro.14O4, usually referred to 
as a strongly correlated system, x"(Q,'^) at incommen- 
surate wave vectors exhibits almost the same magnitude 



and a similar w-dependence |22 . 

In Sr2Ru04, electronic correlations are incorporated in 
RPA calculations: the spin susceptibility xil, ^) becomes 
enhanced through the Stoner-factor I{q) (|,|l^: 
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FIG. 4. Results from fits to a Gaussian profile of 6.2 meV 
constant-Lj scans at Qo=(1.3,0.3,0) along the (0,1,0): temper- 
ature dependences of (a) x" (Qo, 6.2meV) and (b) the intrinsic 
q-width of the magnetic signal, Aq (FWHM). (c) Comparison 
between ^''{1/TiT) observed by ^''O NMR by Imai et al Q 
(□) and the incommensurate contribution calculated from our 
INS measurements (•). Assuming A=33 kOe//iB |^ the two 
scales in this figure are identical. 
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The q-dependence of the Stoner factor, for an individual 
Ru02 plane, reflects the fact that FM interactions are 
favored over antiferromagnetic interactions in Sr2Ru04 : 
in our units, I{q) = 0.86/(1 -f-0.8(a/7r)2g2) eV (q in A^^) 
[^,|o| . INS results point towards a strong enhancement of 
the spin susceptibility by the Stoner factor (see Eq. (||)), 
such that the system should be close to a magnetic in- 
stability at qo. With x'(9o,0) — 180 ^^.eV^\ one de- 
duces from Eq. (|) that 2(^X0(90,0) ~0.99, instead 
of being larger than 1 for a magnetic instability. Thus, 
incommensurate spin fluctuations are stronger than FM 
fluctuations in Sr2Ru04, as suggested in ref. jlOt . 

The temperature dependence of both x" (Qo, 6.2meV) 
and the intrinsic q-width are reported in Fig. ^, as 
deduced from constant-w scans performed at 6.2 meV 
around Q=(1.3,0.3,0) along the (0,1,0) direction at dif- 
ferent temperatures, x" (Qo, 6.2meV) exhibits a sharp 
decrease upon temperature increase and simultaneously 
the magnetic response broadens (the width of the signal 
can be reliably determined only up to 200 K). The T- 
dependence of x" observed in INS measurements may be 
described by the out-smearing of the Fermi surface due 
to thermal hopping of electrons into unoccupied states 
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(see the numerator in Eq. (||)), yielding a lowering of the 
dynamical susceptibility at qo and its broadening in q- 
space. The T-dependence of the magnetic response at qo 
can indeed be qualitatively reproduced psf using Eqs. 

^-iid a description of the LDA band structure by 
three mutually non- hybridizing tight-binding bands 

INS measurements point out the existence of strong 
magnetic response at qo, but do not reveal any sizeable 
FM fluctuations. In contrast, the uniform spin suscep- 
tibility and the Knight shift measurements [ p^j24| 
provide evidence of strong FM correlation in Sr2Ru04. 
However, the delicate balance between FM and incom- 
mensurate spin fluctuations should become visible in the 
spin-lattice relaxation rate Ti measured by both -'^^O and 
i°iRu NMR experiments These NMR-techniques 
probe the low energy spin fluctuations (cj with re- 
spect to INS measurements); furthermore, they integrate 
the fluctuations in q-space. Since the INS studies have 
determined the incommensurate fluctuations on an abso- 
lute scale we may estimate their contribution to (l/TiT), 
^^^(l/TiT). In general (l/TiT) probes the q-summation 
of the the imaginary part of the susceptibility divided by 
the frequency in the limit lo —>■ (i.e. its initial slope), 
^ L~>o; its temperature dependency is shown in 
Fig. 4.C (left scale). 

What renders the quantitative comparison between the 
INS and NMR-results more difficult is the estimate of 
the hyperfine field whose q-dependent Fourier transform, 
A{q), weights the susceptibihty in NMR-studies. Consid- 
ering that INS magnetic fluctuations are sharply peaked 
around qo, one may approximate A(q)~A{qQ), and gets 
0, 
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with \A{qo)\'^ = Afjl + l/2(cos(27r0.3) -I- cos(27r/a0.3))] 
(A = 33 kOe//iB ||l) for and A{qo)=-299 kOe/fiB 
|pl] , H for ^°^Ru. Using these values, we directly com- 
pare ^^■5(l/TiT) with the measured ^^O (l/TiT) in 
Fig. jijc (right scale). Clearly, the spin fluctuations 
at qo significantly contribute to ^^(1/TiT), and can ex- 
plain a large part of the reported T-dependence |]l^ , p4| . 
Similar calculation for the ^'^^Ru (1/TiT) (not shown) 
yields even a stronger contribution. The remaining parts 
in i'^'i°i(I/riT), likely associated with FM excitations, 
should exhibit a less pronounced T-dependence similar 
to that of the uniform static spin susceptibility. Fur- 
thermore, assuming a weak q-dependence for these FM 
excitations the comparison of NMR and INS mea- 
surements allows us to estimate the ferromagnetic charac- 
teristic energy to be of the order of 50 meV. This rather 
elevated value actually provides a satisfactory explana- 
tion for the absence of FM fluctuations in our INS mea- 
surements. 



To conclude, our INS measurements demonstrate the 
existence of incommensurate spin fluctuations related to 
dynamical nesting properties of the Sr2Ru04 Fermi sur- 
face. Our data suggest that the system is close to a 
magnetic instability at qo// = (±0.67r/a,±0.67r/a). The 
comparison of INS and ^''(l/TiT) measurements suggests 
that the FM fluctuations are transferred to higher energy 
with respect to the spin fluctuations at qo. All these re- 
sults cast some doubt on the predominant role of FM 
spin fluctuations in the mechanism of superconductivity 
in Sr2Ru04. 

We wish to acknowledge P. Pfeuty, J. Bobroff and Ph. 
Mendels for helpful discussions. 



[1 
[2 

[s: 
[4: 

[5 
[G 
[7 

[» 

[9 

[lo; 

[11 
[12 

[13 
[14 

[15 
[16 
[17 
[18 

[19] 



[20] 



[21] 



[22] 
[23] 
[24] 
[25] 



Y. Maeno et al, Nature 372, 532 (1994). 
Y. Maeno et al.,J. Phys. Soc. Jpn. 66, 1405 (1997). 
J. J. Neumeier et al, Phys. Rev. B 50, 17910 (1994). 
T. M. Rice et al, J. Phys. Cond. Matt. 7, L643 (1995). 
I. I. Mazin et al, Phys. Rev. B 56, 2556 (1997). 
K. Machida et al, J. Phys. Soc. Jpn. 65, 3720 (1996). 
T. Oguchi, Phys. Rev. B 51, 1385 (1995). 
I. I. Mazin et al, Phys. Rev. Lett 79, 733 (1997). 
A. P. Mackenzie et al, Phys. Rev. Lett. 76, 3786 (1996). 
I. I. Mazin and D.J. Singh, Cond-mat/9902193. 
K. Ishida et al, Phys. Rev. B 56, R505 (1997). 
M. Luke et al, Nature 394, 558 (1998). 
Ishida et al, Nature 396, 658 (1998). 
Imai et al, Phys. Rev. Lett. 81, 3006 (1998). 
Braden et al, Phys. Rev. B 58, 847 (1998). 
Nakastuji et al, J. Phys. Soc. Jpn. 66, 1868 (1997). 
Braden et al, Phys. Rev. B 57, 1236 (1998). 
International Tables of Crystallography, Vol. C, eds A. J. 
C. Wilson, Kluwer Academic publisher (1995). 
S. W. Lovesey, Theory of Neutron Scattering from Con- 
densed matter. (Clarendon, Oxford, 1984), Vol.2, the 
kinematic factor kp/ki has been omitted in Eq. (1) for 
sake of simplicity. 

The spin susceptibility is associated with fluctuations 
of a single spin component: x(Qi'^) — X"'' — 
-{gtiB?^ dte^p{iut){[Sf(t),S^^]), where a,l3 are 
Cartesian coordinates, x is half of the transverse spin 
susceptibility • 

We have used acoustic phonons in scans (2+^,0,0) and 
(2,^,0) for hjjj= 6.75 and 8.26 meV (sound velocity ~ 
20.8 and 43.9 meV.A, respectively) jlj. 

G. Aeppli et al. Science, 278, 1432 (1997). 
P. Pfeuty, private communication. 

H. Mukuda et al, J. Phys. Soc. Jpn. 67, 3945 (1998). 
C. Berthier et al, J. Phys. I France 6, 2205 (1996); R.E. 
Walstedt et al, Phys. Rev. Lett. 72, 3610 (1994). 

[26] For a magnetic field applied along the c axis, A^=|^// 

with A//=-18.5 kOe//ifl and A^=-2Q.9, kOe/iiB 



4 



